PROOF OF THE REFINED ALTERNATING SIGN MATRIX CONJECTURE 



Boron ZEILBER GER 1 

Abstract. Mills, Robbins, and Rumsey conjectured, and Zeilberger proved, that the number of 
alternating sign matrices of order n equals A(n*j : — n\(n-\-\)\ ^ (2n—i)\ * Mills, Robbins, and Rumsey 
also made the stronger conjecture that the number of such matrices whose (unique) '1' of the 
first row is at the r th column, equals A(n) ( n ^Y 2 ) ^n-i^ / ^n-i) ■ Standing on the shoulders of 
• A.G. Izergin, V. E. Korepin, and G. Kuperberg, and using in addition orthogonal polynomials and 

g-calculus, this stronger conjecture is proved. 

INTRODUCTION 



An alternating sign matrix, or ASM, is a matrix of O's, l's, and — l's such that the non-zero elements 
in each row and each column alternate between 1 and — 1 and begin and end with 1, for example: 



o 

U 

C3 



Mills, Robbins, and Rumsey [MRR1][MRR2]([S], conj. 1) conjectured, and I proved[Z], that there 
are 

(N ■ . . l!4!7!---(3ra-2)! 




^3 ■ v ' n\(n + 1)\ ■ ■ ■ (2n - 1)! 

alternating sign matrices of order n. Another, shorter, proof was later given by Greg Kuperberg[K]. 
Kuperberg deduced the straight enumeration of ASMs from their weighted enumeration by Izergin 
and Korepin[KBI]. In this paper, I extend Kuperberg's method of proof to prove the more general, 
refined enumeration, also conjectured in [MRR1] [MRR2] , and listed by Richard Stanley [S] as the 

S' 
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third of of his "Baker's Dozen", that: 
Main Theorem.— There are 



Mn,r) := A(n) K } 

\ n-l ) 

n x n alternating sign matrices for which the (unique) '1' of the first row is at the r th column. 

As in Kuperberg's proof, we are reduced to evaluating a certain determinant. Unlike the original 
determinant, it is not evaluable in closed form. We evaluate it by using the q-analog of the Legendre 
polynomials over an interval, introduced, and greatly generalized, by Askey and Andrews[AA] and 
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Askey and Wilson [AW]. All that is needed from the general theory of orthogonal polynomials and 
from g-calculus is reviewed, so a sufficient condition for following the present paper is having read 
Kuperberg's paper [K] that includes a very clear exposition of the Izergin-Korepin formula, and of 
its proof. Of course, this is also a necessary condition. In order to encourage readers to look up 
and read Kuperberg's beautiful paper [K], and to save myself some typing, I will use the notation, 
and results, of [K], without reviewing them. 2 

Boiling It Down To a Determinant Identity 

Let B(n,r) be the number of ASMs of order n whose sole '1' of the first (equivalently last) row, is 
at the r th column. In order to stand on Kuperberg's shoulders more comfortably, we will consider 
the last row rather than the first row. 

Like in [K], define [x] := (q x / 2 — q~ x / 2 ) / (q 1 / 2 — g -1 / 2 ), take q := e 2 v-i 7r / 3 j anc ; consider 

Z(n; 2,..., 2, 2 + a; 0,...,0) , 

where between the two semi-colons inside Z there are n-12's followed by a single 2 + a. Here a 
is an indeterminate. Let's look at an ASM of order n, whose sole '1' of the last row is at the r th 
column. It is readily seen that the r — 1 zeros to the left of that '1' each contribute a weight of 
[2 + a], while the remaining n — r zeros, to the right of the aforementioned '1', each contribute a 
weight of [1 + a]. The '1' itself contributes q~ l ~ a / 2 , which is q~ a / 2 times what it did before. Hence 

n 

Z(n; 2,..., 2, 2 + a; 0,...,0) = (-l) n q~ n q~ a/2 ^ B(n,r)[2 + a] r_1 [l + a] n ~ r . 

r=l 

We also have, thanks to [K], (or [Z]:just plug in a = above): 

Z(n;2,...,2,2;0,...,0) = (-l)»<TM(n) . 

Hence: 
Since 

{[2 + a] r - 1 [l+a} n - r ; l<r<n} 

are linearly independent, the B(n,r) are uniquely determined by the above equation. Hence the 
Main Theorem is equivalent to: 

Z(n;2,...,2,2 + fl;0,...,0) _ q~ a / 2 A (n + r - 2\ (2n - 1 - r\ x 
Z(n;2,...,2,2;0,...,0) ~ ^ n - 1 ) { n-1 ) [2 + ^ [1 + ° J " 



2 

Kuperberg's paper may be obtained from his Home Page ittp : //www .math. yale . edu/users/greg. Once you down- 
loaded the file asm.ps.gz you should first get rid of the .gz by typing gunzip asm.ps, and then to get a copy, you 
type lpr asm.ps. 
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By replacing re by re + 1, and changing the summation on r to start at 0, we get that it suffices to 
prove: 



Z(n + l;2,...,2,2 + q;0,...,0) _ g~ a / 2 ^ (n + r\ (2n - r\ 
Z(re + 1;2,...,2,2;0,...,0) " («) ^ V n A - + aj[l + «j 



Let Z(n; x±, . . . , x n ; yi, . . . , y n ) denote the right hand side of the Izergin-Korepin formula (Theorem 
6 of [K]). First replace n by re+1. Then, taking Xj = 2 + ie, for i = 1, . . . , re, x n+ i = 2 + a + (re+ l)e, 
and j/j = — (j — l)e, j = 1, . . . , re + 1, yields, after cancellation, 

Z(n + 1 ; 2 + e, . . . , 2 + ree , 2 + a + (re + l)e ; 0, -e, . . . , -ree) 
Z(n + 1 ; 2 + e, . . . , 2 + ree , 2 + (re + l)e ; 0, -e, . . . , -ree) 



-a/2 -rj [2 + a + (re + 1 + j)e] [1 + a + (re + 1 + j)e] Vr [(re - j)e] detM ra+1 (a) 
9 1 = 1 [2 + (re + 1 + j)e] [1 + (re + 1 + j)e] " ' f} Q [a + (re - j)e] ' detM n+1 (0) 

where M re+ i(a) = {rriij , < i, j < re) is the (re + 1) x (re + 1) matrix, defined as follows: 



_ J l/([2 + (i+j + l)e][l + (i + j + l)e]), if < i < re - 1, < j < re; 

l/([2 + a + (re + J + l)e][l + a+ (re + j + l)e]) , if i = re, < j < re. 



Taking the limit e — > 0, replacing g e by s, g a by X, setting w := e v/ ~ T7r / 3 , evaluating the limit 
whenever possible, and cancelling out whenever possible, reduces our task to proving the following 
identity: 



lim 



(l-s)"detAT n+1 pQ | 
i\ detiV n+ i(l) J 



2 \r\n — r 



-(l-X) n {^3) n+2 w- n ^ _ r (n + r\(2n-r\. 2 
n!(l + X + X2)«+i( 3 ' l + 1 ) ^ V n )\ re ^ ;i > 

(NotYetDone) 

Here the matrix N n+1 (X) is M n+ i(a) divided by 3, to wit: N n+1 (X) = (pij,0 < i,j < re) is the 
(re + 1) x (re + 1) matrix, defined as follows. For the first re rows we have: 

Pi * = i _ S 3(i+j+i) ' (0<i<n-l , 0<j<re) , 
while for the last row we have: 

1 - Xs n+j+1 

Pn > j= l-X3 8 3(n+j+l) ' • 

We are left with the task of computing the determinant of N n+1 (X), or at least the limit on the 
left of (NotYetDone). 
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A Short Course on Orthogonal Polynomials 



I will only cover what we need here. Of the many available accounts, Ch. 2 of [Wilf] is especially 
recommended. For the present purposes, section IV of [D] is most pertinent. 

Theorem OP.— Let T be any linear functional ('umbra') on the set of polynomials, 
and let c, := T(x l ) be its so-called moments. Let 



A n := det 



/ c Cl 
Cl c 2 



V 



Cn+1 



C2r. 



If A n ^ 0, for n > 0, then there is a unique sequence of monic polynomials P n (x), where the degree 
of P n (x) is n, that are orthogonal with respect to the functional T: 

T(P n (x)P m (x))=0 if m + n . 

Furthermore, these polynomials P n (x) are given 'explicitly' by: 



Pn{x) 



1 



/ C Ci 

Ci c 2 



A 



det 



n-l 



Cn — 1 Cjj 



V i 



X X 



Cn+l 



Proof: □. 



Corollary 1.— 



Proof: □. 



T{x n P n {x)) 



for n > 1 



Corollary 2.— If S 1 is another linear functional, <ij := S{x l ) , and 



r„ := det 



/c Cl ... 

Cl c 2 ... 

Cn — 1 C n . . . 

\ d di d 2 



Cn+1 

C2n-1 
d n ' 



Then 



Proof: □. 



r n _ S{P n {x)) 
A n T{x n P n {x)) 



for n > 1 



(GeneralFormula) 



For a long time it was believed that theorem OP was of only theoretical interest, and that, given the 
moments, it was impractical to actually find the polynomials P n (x), by evaluating the determinant. 
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This conventional wisdom was conveyed by Richard Askey, back in the late seventies, to Jim 
Wilson, who was then studying under him. Luckily, Wilson did not take this advice. Using 
Theorem OP lead him to beautiful results [Wils], which later lead to the celebrated Askey- Wilson 
polynomials [AW]. Jim Wilson's independence was later wholeheartedly endorsed by Dick Askey, 
who said: "If an authority in the field tells you that a certain approach is worth trying, listen to 
them. If they tell you that a certain approach is not worth trying, don't listen to them". 

The 'uselessness' of (GeneralFormula) was still proclaimed a few years later, by yet another 
authority, Jean Dieudonne, who said ([D], p. 11): "La formule generale [(GeneralFormula)] 
donnant les . . . sont impraticables pour le calcul explicite . . . ". 

There is another way, by which (GeneralFormula) , and its immediate corollaries 1 and 2 could 
be useful. Suppose that we know, by other means, that a certain set of explicitly given monic 
orthogonal polynomials Q n (x) are orthogonal with respect to the functional T, i.e. T(Q n Q m ) = 
whenever n / m. Then by uniqueness Q n = P n . If we are also able, using the explicit expression 
for Q n (x), to find T(x n Q n (x)), then Corollary 1 gives a way to explicitly evaluate the Hankel 
determinant A n . If we are also able to explicitly compute S(Q n (x)), then we would be able to 
evaluate the determinant T n . This would be our strategy in the evaluation of the determinants on 
the left of (NotYetDone), but we first need to digress again. 

A Lean and Lively Course in q-Calculus 

Until further notice, 

(a) n := (1 - o)(l - qa)(l - q 2 a) ... (1 - q^a) . 



If I had my way, I would ban 1— Calculus from the Freshman curriculum, and replace it by 
q— Calculus. Not only is it more fun, it also describes nature more accurately. The traditional 
calculus is based on the fictitious notion of the real line. It is now known that the universe is quan- 
tized, and if you are at point x, then the points that you can reach are in geometric progression 
q % x, in accordance with Hubble expansion. The true value of q is almost, but not quite 1, and is a 
universal constant, yet to be determined. 



The g-derivative, D q , is defined by 

D q f(x) 

The reader should verify that 



fix) - f(qx) 
(1 - q)x 



1 — a a 
D q x a = — — i"" 1 



1-q 



and the product rule: 

D q [f(x) ■ g(x)\ = f(x) ■ D q g(x) + D q f(x) ■ g(qx) . (ProductRule) 
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The q-analog of integration, independently discovered by J. Thomae and the Rev. F.H. Jackson 
(see [AA], [GR]), is given by 



/ f(x)d q x := a(l - q)y2f(aq r )q r , 
Jo r=0 



and over a general interval: 



pd pd pc 

/ f(x)d q x:= / f(x)d q x- / f(x)d q x 

Jc JO Jo 



The reader is invited to use telescoping to prove the Fundamental Theorem of q-Calculus 

cd 



D q F(x) d q x = F{d) - F(c) . 



Combining the Product Rule and the Fundamental Theorem, we have: 
q-Integration by Parts: If f(x) or g(x) vanish at the endpoints c and d, then 

i-d nd 

/ f(x) ■ D q g(x)d q x = - / D q f{x) ■ g{qx) d q x . 

J c J c 

Corollary: If g(q l x) vanish at the endpoints c and d, for i = 0, 1, . . . , n — 1, then 
[ f{x)-D n q g{x)d q x = {-ir f D n q f{x)-g{q n x)d q x . 

J c J c 

Even those who still believe in 1-Calculus can use g-Calculus to advantage. All they have to do is 
let q — > 1 at the end. 

The q-analog of J s x a dx = (1 — s a+1 )/(a + 1) is 

If 1 l — s a + 1 

/ x a d a x = — . (q — Moment) 

l-qj s q 1 

Now this looks familiar! Letting q := s 3 in the definition of N n+1 (X), given right after (NotYetDone), 
(this new l q' has nothing to do with the former Kuperberg q), we see that the matrix N n+ i(l) is 
the Hankel matrix of the moments with respect to the functional 

i r 1 



T(f(x)) := -L_ f f( x ) dqX 



So all we need is to come up with orthogonal polynomials with respect to the 'q-Lebesgue-measure', 
over the interval [s,l]. 
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q-Legendre Polynomials 

The ordinary Legendre polynomials, over an interval (a, b) may be defined in terms of the Rodrigues 
formula 

—D"{(x-ar(x-br} . 
The orthogonality follows immediately by integration by parts. This leads naturally to the g-analog, 

Q n (x; a, b) := ±—^-D^{(x - a)(x -qa)...(x- aq^ 1 ) ■ (x - b)(x - qb) . . . (x - bq^ 1 )} . 

Using ^-integration by parts repeatedly, it follows immediately that the Q n (x;a,b) are orthogonal 
w.r.t. to q— integration over (a, b). The classical case a = —1, 6 = 1 goes back to Markov. Askey 
and AndrewsfAA] generalized these to q-Jacobi polynomials, and Askey and WilsonfAW] found the 
ultimate generalization. While at present I don't see how to apply these more general polynomials 
to combinatorial enumeration, I am sure that such a use will be found in the future, and all 
enumerators are urged to read [AA], [AW], and the modern classic [GR]. 

Going back to the determinant N n+1 (X) of (NotYetDone), we also need to introduce the functional, 
defined on monomials by: 

1 _ Xs n+j+1 
S ( X ^ = 1 - x 3 q n +3+ 1 ' 

and extended linearly. 

Let X := q a / 3 . Then (recall that s = q 1 / 3 ): 

S(x J ) = ^-^r = / x a+n x j d q x . 

y ' 1 - qn+j+a+i 1-q J s q 

By linearity, for any polynomial p(x): 

i r 1 

S(p(x)) = —— / x a+n p(x) d q x . 



Using Corollary 2 of (General Formula) , we get 

det N n+1 {X) x a+n P n (x) d q x 



(AlmostDone) 



det7V n+1 (l) " f\ n P n (x)d q x 
where P n {x) is now the q— Legendre polynomial over [s, 1], Q n (x; s, 1) and s = q 1 ^ 3 . In other words: 

Pn(x) := fe|^^ n {(x - l)(x -q)...(x- q^ 1 ) ■ (x - s)(x - qs) . . . (x - sq^ 1 )} . 



(q n+1 ) 



n 
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DENOUEMENT 



It remains to compute the right side of (Almost Done). Let's first do the denominator. 
Proposition Bottom.— 



l- q J s X ^ n[X) aqX (qn+i 



a (q) 2 n (q- n s) 2n+1 
) n (q n + 1 ) n+1 

First Proof: Use g-integration by parts, n times (i.e. use the above corollary). The resulting 
g-integral is the famous g-Vandermonde-Chu sum, that evaluates to the right side. See [GR], or 
use qEKHAD accompanying [PWZ]. □ 

Remark: Proposition Bottom, combined with Corollary 1 of (GeneralF ormula) gives an alterna- 
tive evaluation of Kuperberg's determinant N n (l), needed in [K]. 

Second Proof: Don't get off the shoulders of Greg Kuperberg yet. Use his evaluation, and 
Corollary 1 of (GeneralF 'ormula). .□ 

Proposition Top.— Recalling that X = q a ^ 3 , we have 

-L- f x^P n (x) d q x = ^J?* 3 ^ ■ f; jf ( 9 «+* - f)(q^ k - f+W) 



( l) n (qX 3 ) n ^ ^2^ k+1/3x3k+1 -j-J^n+fc+l/3 _ ^^n+k+1/3 _ q r+l/3} 



ln k=0 r=0 



Proof: Let 



fl — aY 

Fn(x) := t^TlH* ~l)(x-q)...(x- q n ~ l ) ■ (x - s)(x - qs) . . . (x - sq n ~ l ) , 

\q )n 

so that P n (x) = D™F n (x). Since F n (q l x) vanish for i = 0, . . . , n — 1 at both x = 1 and x = s, we 
have by g-integrating by parts n times (the above corollary), that 

/ x n+a ■ Pn(x) d q x = ! x n+a ■ D n q F n (x) d q x = (-l) n f D n q {x n+a ] ■ F n (q n x) d q x 

J s J s J s 

= t jr s ^f- a ^( q n x)d q x . 

Now use the definition of ^-integration over [s, 1] and replace q a by X 3 , to complete the proof. □ 

To compute the right side of (AlmostDone), we only need to divide the expression given by 
Proposition Top by the expression given by Proposition Bottom. Doing this, multiplying by 
(1 — s) n = (1 - q l / 3 ) n , and taking the limit q — > 1, we get that the left side of (NotYetDone) 
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is (Warning, now we are safely back in 1— land, so from now (a) n := a(a + 1) . . . (a + n — 1), the 
ordinary rising factorial): 



After trivial cancellations, equation (NotYetDone) boils down to 



(-l)"(3n + l)! ^ 
( 1-A } 3 ^n!3(-n + l/3) 2n+1 2J* + lM* + 2/^»* 



= (\^) n E »" r_ " ^ r ) ~ r j (1 + wl) r (1 - w 2 X)"" r . (Done) 

This was given to EKHAD, the Maple package accompanying [PWZ]. EKHAD found a certain linear 
homogeneous second order recurrence in n that is satisfied by both sums on the left of (Done) (and 
hence by their difference), and also by the right side. It remains to prove that both sides of (Done) 
agree at n = 0, 1, which Maple did as well, even though it could be done by any human. 

The input and output files are obtainable by anonymous ftp to ftp.math.temple.edu (login as 
anonymous), cd-ing to directory /pub/zeilberg/ref ined, and geting files inDone and outDone 
respectively. EKHAD (and qEKHAD) are available at directory /pub/zeilberg/programs. Alterna- 
tively, on the Web, go to my Home Page given at footnote 1, and click in the appropriate places. 
Of course, your computer should be able to reproduce file outDone. Once you have downloaded 
EKHAD and inDone into a directory, type: maple -q<inDone>outDone. After 380 seconds of CPU 
time, outDone would be ready. 
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